Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN WF$t style file v2.2) 



On gravito-magnetic time delay by extended lenses 
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ABSTRACT 

Gravitational lensing by rotating extended deflectors is discussed. Due to spin, cor- 
rections on image positions, caustics and critical curves can be significant. In order 
to obtain realistic quantitative estimates, the lens is modeled as a singular isother- 
mal sphere. Gravito-magnetic time delays of ~ 0.2 days between different images of 
background sources can occur. 
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1 INTRODUCTION 

Mass-energy currents relative to other masses generate 
space-time curvature. This phenomenon, known as intrin- 
sic gravito-magnetism, is a new feature of general theory 
of relativity and other conceivable post-Ne wtonian theo- 
ries of gravity Csee lCinfolini fc Wheelerlll 99,4 and references 
therein). The gravito-magnetic field has not been yet de- 
tected with high accuracy. Some results have been reported 
from laser-ranged satellites, when the Lense-Thirring pre- 
cession, due to the Earth spin, was measured by studying 
the o rbital perturbations of LAGEO S and LAGEOS II satel- 
lites ijCiufolini fc Pavlisl ll998. 2004J)- According to a prelim- 
inary analysis, the predictions of the general theory of rela- 
tivity have been found to agree wit h the experimental value s 
within the ~ 10 per cent accuracy llCiufolini fc Pavlii2004l) . 
The NASA's Gravity Probe B satellite should improve this 
measurement to an accuracy of 1 per cent. Intrinsic gravito- 
magnetism might play a relevant role also in the dynam- 
ics of the accretion disk of a supermassive black hole or in 
the alignment of jets in ac tive galactic nuclei and quasars 
jCiufolini fc Wheelerlll99l . 

Whereas the tests just mentioned limit to the gravi- 
tational field outside a spinning body, the general theory of 
relativit y predicts peculi ar phenomena also inside a rotating 
shell fsee lWeinbersll972l for example). Gravitational lensing 
can represent a to ol to fu l ly test the effects of the gravito- 
magnetic field fsee lSerenol (l2003af) and reference therein). In 
this paper, we are interested in the gravito-magnetic time 
delay induced in different images of the same source due to 
gravitational lensing. Whereas here we are mainly faced with 
intrinsic gravito-magnetism, i.e. with spinning deflectors, a 
translational motion of the lens can also induce interesting 
phenomena, which in the framework of general relativity are 
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strictly connected to the Lorentz transformation properties 
of the gravitational field. The effect of the deflector's velocity 
has been recently observed in the Jovian deflection exper- 
iment conducted at VLBI, which me asured the time delay 
of light from a background quasar iFomalont fc Kopeikinl 
2003). Although a controversy has emerge d over the theo - 
retical interpretation of this measurement <lKopeikinll2004t) . 
there is agreement about the role of the deflector's motion. 

Gravitational time delay by spinning deflectors has 
been addressed by several autho rs with very differ- 
ent approaches. iDvmnikoval lll986l) discussed the addi- 
tional time-delay due to rotation by integrating the light 
geodesi es of the Kerr m etric. Using the Lense-Thirring 
metric, iGlicensteirJ lll999l) considered the time delay for 
light rays pas sing outside a s p inning star. Kopeikin and 
collaborators jKopeikinl Il997t iKopeikin fc Schaferl ll99St 
iKopeikin fc Mashhoo J2 002 ) analysed the gravito-magnetic 
effects in the propagation of light in the field of self- 
gravitating spinning bodies. The gravitational time de- 
lay due to r o tating masses was fur t her d iscussed in 
ICiufolini etafl feOOST ; ICiufolini fc Riccil (120031) . where the 
cases of light rays crossing a slowly rotating shell or propa- 
gating in the field of a distant source were analyzed in the 
linear approximation of general theory of relativity. Effects 
of an intrinsic gravito-magnetic field were further st udied in 
the usual framework of gravitational le nsing theory iSerenol 
l2002ll2003b1lal:ISereno fc CardonefeOO^ . i.e. i) weak field and 
slow mot ion approximation for the lens and it) thi n lens hy- 
pothesis (Schneider ct al. 19921: iPetters et al .120011) . Expres- 
sions for bending and time delay of electromagnetic waves 
were found for sta tionary spinn ing deflectors with general 
mass distributions <ISerenoll2002r) . 

The paper is as follows. In Section basics of gravi- 
tational lensing by a stationary extended deflector are re- 
viewed. In Section [3 we introduce our reference model for 
the lens, i.e. a spinning singular isothermal sphere. Rele- 
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vant lensing quantities for such a deflecting system are de- 
rived in Section [I] In Section |KJ the lens equation is solved. 
A quantitative discussion of the gravito-magnetic time de- 
lay is in Section |S| where we also investigate the effect on 
the determination of the Hubble constant. Section |7| is de- 
voted to some final considerations. Unless otherwise stated, 
throughout this paper we consider a flat cosmological model 
of universe with cosmological constant and a pressureless 
cosmological density parameter Qmo = 0.3 a Hubble con- 
stant H = 72 km s^Mpc" 1 . 



2 GRAVITO-MAGNETIC DEFLECTION 
POTENTIAL 

Gravitational lensing theory can be easily developed in 
the gravitational field of a rotating stationary sourc e when 
the linear approximation of general relativity holds (ISerenol 
2002). The time delay of a kinematically possible ray, with 
impact parameter £ in the lens pla ne, relative t o the un- 
lensed one is, for a single lens plane . (ISerenoll2002T) 
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where ?/> is the deflection potential; D s , Dd and Dds are 
the angular diameter distances between observer and source, 
observer and lens and lens and source, respectively; Zd is the 
lens redshift; rj is the bidimensional vector position of the 
source in the source plane. We have neglected a constant 
term in equation ffl . since it has no physical significance 
JSchneider et alJll992l) . 

The deflection potential can be expressed as the sum of 
two terms 



ip~tpo + V'grm; 
the main contribution is 
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where £o is a length scale in the lens plane and E is the 
projected surface mass density of the deflector, 
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the gravito-magnetic correction to the d eflection pote ntial, 
up to the order v/c, can be expressed as (Scrcno 2002) 
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where (v-ei n ); is the weighted average, along the line of sight 
ei n , of the component of the velocity v along ei n , 



JWf») Pit I) dl 

£(€) 



(6) 



In the thin lens approximation, the only components of the 
velocities parallel to the line of sight enter the equations of 
gravitational lensing. We remind that the time delay func- 
tion is not an observable, but the time delay between two 
actual rays can be measured. Similar results, based on a 
multi-polar description of t he gravitational field of a sta- 
tionary lens, can be found in iKopeiki 



3 SINGULAR ISOTHERMAL SPHERE 

Isothermal spheres are widely used in astrophysics to under- 
stand many properties of systems on very different scales, 
from galaxy haloes t o clusters of galaxies iMo et al J 11993 
ISchneider et al]ll992T) . In particular, on the scale relevant to 
interpreting time delays, isot hermal models are favoured b y 
data on early-type galaxies jKochanek fc Schechte'ill2004h . 
The density profile of a singular isothermal sphere (SIS) is 



p(r) = 



(7) 



2ixGr 2 ' 

where <r v is the velocity dispersion. The corresponding pro- 
jected mass density is 

E«) = if (8) 

Since the total mass is divergent, we introduce a cut-off 
radius Rsis 2> £• The cut-off radius of the halo must be much 
larger than the relevant length scale which characterizes the 
phenomenon, in order to not significantly affect the lensing 
behavior. The total mass of a truncated SIS is 

2cr 2 

M SIS = -£-R S is- (9) 

A limiting radius can be defined as r n , the radius within 
which the mean mass density is n times the critical density 
of the universe at the redshift of the galaxy, Zd- For a SIS, 
it is 

2(j v 
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where H is the time dependent Hubble parameter. 

The total angular momentum of an halo, J, can be ex- 
pressed in terms of a dimensionless spin parameter A, which 
represents the ratio between the actual angular velocity of 
the system and the hypothet ical angular velocity that is 
needed to support the system jPadmanabhanir2002l) . 



J = A 
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where M and E are the total mass and the total energy 
of the halo, respectively. In the hypothesis of initial an- 
gular momentum acquired from tidal torquing, typical val- 
ues of A can be obtained from the relation between energy 
and virial radius and the details of spherical top-hat model 
ljPadmanabhanl2002l 'l . As it was derived from numerical sim- 
ulations, the distribution of A is nearly independent of the 
mass and the power spe ctrum. It can be appr oximated by a 
log-normal distribution l|Vitvitska et al]l2002f) 
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with A ~ 0.05 and o\ ~ 0.5. The distribution peaks around 
A ~ 0.04 and has a width of ~ 0.05. 

Fr om the virial th eorem, the total energy is easily ob- 
tained jMo et al.lll993l 



Esm = — Msisfv 



(13) 



Finally, the total angular momentum of a truncated SIS can 
be written as 



Jsis = A- 



G 



(14) 
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In general, the angular velocity of a halo is not con- 
stant and a differenti al rotation should be considered 
l|CaDozziello et aDl2003t) . However, assuming a detailed ro- 
tation pattern does not affect significantly the results. In 
what follows, we will consider the case of constant angular 
velocity. 



4 LENSING BY A ROTATING SIS 

Let us consider gravitational lensing by a SIS in rigid rota- 
tion, i.e with a constant angular velocity a;, about an arbi- 
trary axis passing through its center. The deflection angle 
can be written as (|Sereno fc Cardonel l2002'l 
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where (£, i9) are polar coordinates in the lens plane and Wi 
and W2 are the components of uj along the fi- and the £2- 
axis, respectively. The parameter u has to be interpreted 
as an effective angular velocity uj = Jsis//sis where 7sis is 
the central momentum of inertia of a truncated SIS, /sis = 
2/9Afsis-Rsis- in terms of the spin parameter, 



u = 9A 



Rsi 
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E quations l|15J and ILLSI correct the result in equation (9) 
in ICapozziello fc Ra (|2001j), which was obtained under the 
same assumptions but was affected by an error in the compu- 
tation of the integrals. For a SIS, there are two main contri- 
butio ns to the gravito-magneti c correction to the deflection 
angle iSereno fc Cardonel2002r) : the first contribution comes 
from the projected momentum of inertia inside the radius 
£; the second contribution is due to the mass outside £ and 
can become significant in the case of a very extended lens, 
i.e. for a very large cut-off radius. 

Let us consider a sphere rotating about the ^2-axis, 
uj\ — 0,^2 = u. In order to change to dimensionless vari- 
ables, we introduce a length scale, 

2 D d D ds 
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The dimensionless position vector in the lens plane is x 
£/£o- The dimensionless deflection potential ip, 

D d D ds 
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can be written as 
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where x = |x| and L = (2/3)(lL>Re/c) is an estimate of the 
rotational velocity and r is the cut-off radius in units of Re- 
When L > 0, the angular momentum of the lens is positively 
oriented along X2- The dimensionless Fermat potential <f> is 
defined as, 



x ,y) = o( x ~y) 
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where y = n/ (^o\- 

The scaled deflection angle is related to the dimension- 
less gravitational potential ip 



a(x) = V^(x). 
We get 
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The determinant of the Jacobian matrix reads 
AsisixxiXi) ~ 1- — - L— - - V 

X X \ X) 

The convergence k = E/E cr is the projected density in units 
of the critical surface density, 

S cr = C „ . (26) 
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fesis is positive when x\ > —1/3L. Since we have introduced 
a cut-off radius, this condition holds for x\ ~ r. The condi- 
tion Lr < 1/3 warrants that fcsis > for all points in the 
lens plane. 



5 THE LENS EQUATION 

In general, the inversion of the lens equation is a math- 
ematical demanding problem. However, since the gravito- 
magnetic effect in the weak field limit is an higher-order 
correction, a rotating system can be studied in some de- 
tails using a perturbative approach JSerenol2003al) . The lens 
equation can be expressed in a dimensionless form as 



y = x - a(x). 



(28) 



The unperturbed images are solutions of the lens equation 
for L — 0. When y = |y| < 1, a static SIS lens produces 
two images, collinear with the source position and the lens 
centre, at radii x = y + 1 and x — y — 1. When y > 1, only 
one image forms, at x = y + 1. 

Under the condition L <g; 1, we can obtain approximate 
solutions to the first-order in L, given by 

x~x {0 ) +Lx (1) , (29) 

where X( ) and xm denote the zeroth-order solution (i.e. the 
solution of the lens equation for L = 0), and the correction 
to the first-order, respectively. By substituting equation 129> 
in the vectorial lens equation, equation J2HJ, we obtain the 
first-order perturbations, 
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For a source on the j/i-axis (y 2 = 0), equations 13011 and 13H 
reduce to 



U(i)i = + 2 \ x (°n\ 



(32) 
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Figure 1. Image positions (square boxes) and critical line (full 
line) for a rotating SIS lens. The source is the grey circle. The 
centre of the coordinate-axes, the source and the two unperturbed 
images (empty boxes) lie on a straight line. Two images (filled 
boxes) are counter-clockwisely rotated, about the line of sight 
through the centre, with respect to this line. It is L = 10 — 3 and 
r = 50. 
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Figure 2. The shift in the image positions (in units of arcsec) 
due to the gravito-magnetic field for source moving parallelly to 
the j/i-axis. Thick and thin lines are for y 2 = 0.01 and 0.1, re- 
spectively. Full and dashed lines refer to the two images. It is 
L = 4xl0" 5 and r = 15. 
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When L > 0, photons with an impact parameter x± < 
(xx > 0) go around the lens in the same (opposite) sense of 
the deflector. Due to the contribution to the deflection angle 
from the mass outside the impact parameter, photons which 
impact the lens plane on the rri-axis at xi > move closer to 
the centre. This feature, which is peculiar of light rays prop- 
agating inside a halo, is opposite to the case of propagation 
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Figure 3. A source (the grey circle) inside the central caustic 
of a rotating SIS is multiple imaged in a cross shaped pattern; 
the four filled box locate the four images. The critical line is also 
plotted. It is r = 15 and L = 3.5 xlO -5 . 



outside a spinning sphere jSerenol Eo03a). when photons in 
the equatorial plane moving in the same (opposite) sense 
of rotation of a deflector form closer (farther) images with 
respect to the non-rotating case. In fact, in the last case, 
there is no mass outside the photon path. If we do not limit 
to the equatorial plane, for L > 0, the images are rotated 
counter-clockwisely around the line of sight with respect to 
the static case, see Fig.Q 

Let us consider the gravito-magnetic correction to the 
deflection angle for a typical galaxy lens at z& = 0.3 with 
cr v ~ 250 Km s _1 , J?sis ~ 100 Kpc and A ~ 0.1, and a back- 
ground quasar at z s = 2.0. Such a configuration correspond 
to L ~ (2 - 4) x 10" 5 and r ~ 15. For some particular source 
positions, the shift in the image positions in the lens plane 
with respect to the static case can be as large as 0.1". In 
Fig. [5] we plot the shift in the image positions for a source 
moving along the j/i-axis, i.e for 2/2 = const. The maximum 
variation occurs when the source nearly crosses the projected 
rotation axis. 

The critical curve is slightly distorted. The solution of 
det A(xi, xz) — 0, with respect to x%, is 



Xl 



(34) 



The area of the critical curve slightly grows and its centre 
shifts of L along the rri-axis. The critical curve intersects 
the Xi-axis in x\ ~ L±l. The changes in the width and in 
the height of the critical curve are of order 0(L 2 ). 

The four extremal points of the critical curve can be 
mapped onto the source plane through the lens equation 
to locate the cusps of the caustic. We find a diamond- 
shaped caustic with four cusps, centred in {2/1,2/2} = 
{i(§r-l),0}. The axes, of semi- width ~ L 2 , are parallel 
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Figure 4. The gravito-magnetic time delay (in units of days) for 
a source moving parallelly to the j/i-axis. Thick and thin lines are 
for 2/2 = 0.01 and 0.1, respectively. It is L = 4xl0 -5 and r = 15. 
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Figure 5. Relative error in the estimate of the Hubble constant, 
due to neglecting the gravito-magnetic field, for a source moving 
with fixed y 2 = 0.1. It is L = 2.5xl0" 3 and r = 15. 



to the coordinate axes. The orientation and the position on 
the caustic depends on both the strength and orientation of 
the angular momentum and on the radius of the lens. When 
the source is inside the caustic, two additional images form. 
Since the axial symmetry is broken by the gravito-magnetic 
field, the Einstein ring is no more produced. A source, which 
is inside the central caustic, is imaged in a cross pattern, see 
Fig.H 



6 TIME DELAY 

As we have seen in Section |5J the time delay for a spinning 
SIS is made of three main contributions: the geometrical 
time delay, At gcom , the unperturbed gravitational time delay 
by a static SIS, Ato, and, finally, the gravito-magnetic time 
delay, AiGRM, 



At ~ At E 



+ At,, + At GF 



(35) 



The main contribution to the gravitational time delay of a 
light ray with respect to the unperturbed path can be re- 
written as 

2 % (36) 



cAt = 4tt(1 + z d )D d (Jjj 



where 9 is the modulus of the angular position in the lens 
plane, 9 — £/D d . The gravito-magnetic time delay is 



cAtr 



~ -367r(^) (1 + z d )D d 0i\ (38) 

where 9i is the angular distance of an image from the pro- 
jected rotation axis and #sis = Rsis/D d . To obtain equa- 
tion 1381 . we have used the relation i<r. An approximate 
relation holds between Ato and the gravito-magnetic time 
delay, 
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Since L< 1, the angular separation between the two images 
is nearly 



(40) 



and the gravito-magnetic induced retardation between the 
two images, ATgrm = At G RM(x a ) - At G RM(x 6 ), can be 
approximated as 



cATgrm ^ 288tt costf 



X. 



(41) 



Since the two images are nearly collinear with the lens cen- 
tre, the gravito-magnetic time delay is nearly independent 
of the cut-off radius. 

For a typical lensing galaxy at z d = 0.5 with <r v ~ 
250 Km s~ , and a background source at z B — 2.0, the 
gravito-magnetic time delay is ~ 0.1 — 0.2 days for A ~ 0.05- 
0.1. In Fig. 01 the time delay between the images is plotted 
as a function of the source position for a source moving 
perpendicularly to the projected rotation axis. The distance 
of the source from the j/i-axis determines the width of the 
transition between the extremal values. 



6.1 The Hubble constant 

Any gravitational le nsing system c an be used to determine 
the Hubble constant jRefsda]||l964t) . Neglecting the gravito- 
magnetic correction can induce an error in the estimate of 
the Hubble constant. In general, it is 



Ho At = T(a v , 2 d , z s ; Qio), 



(42) 



where the dimensionless function T depends on the lens pa- 
rameters and on the cosmological density parameters, but 
the latter dependence is not very strong. A lens model which 
reproduces the positions and magnifications of the images 
provides an estimate of the scaled time delay H§At be- 
tween the images. Therefore, a measurement of At will yield 
the Hubble constant. Let us consider a rotating galaxy, de- 
scribed by a SIS with known dispersion velocity and redshift, 
which forms multiple images of a background quasar at red- 
shift Zs- An observer can measure the time delay between 
the two images, AT b s , and their positions, x a and X;,. The 
unknown source position y can be obtained by inverting the 
lens equation. In terms of the dimensionless Fermat poten- 
tial d>, the measured Hubble constant turns out to be 
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where 
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and r is the angular diameter distance in units of c/Ho- 

If, we assume a static lens model to model the data, the 
estimated Hubble constant is 
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where the "not correct" estimated position of the source is 
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The relative error in the determination of the Hubble con- 
stant is 

ST _ 

(47) 
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For a source at fixed 3/2, the maximum error is ~ 1/2 \L/y2\, 
see Fig.|5] Since usually L ~ 10 -4 , the induced relative error 
is really negligible. 



7 DISCUSSION 

Since the physics of gravitational lensing is well understood, 
the gravito-magnetic time delay may provide a new observ- 
able for the determination of the total mass and an gular 
momentum of the lensing body jCiufolini fc Riccil2003T) . Al- 
though a detailed model may be required to reproduce the 
overall mass distribution in the lens, interpretation of time 
delay is based on a limited number of parameters. Provided 
the cluster where the deflector lies can be described by a 
simple expansion, the only parameters needed to model the 
time delay are those needed to vary the average surface den- 
sity of the lens near the images and to change the ratio 
between th e quadrupole moment of th e lens and the en- 
vironment jKochanek fc S chechter 2004J). Furthermore, the 
presence of an observed Einstein ring can provide strong 
independent constraints on the mass distribution. 

We have developed our treatment of the gravito- 
magnetic time delay by modeling the lens as a SIS. Isother- 
mal models are supported by both theoretical prejudices 
and estimates from observations of early-type galaxies so 
that SIS turns out to be a surprisingly realistic starting 
point for modeling lens potentials. Gravito-magnetic time 
delays of ~ 0.1-0.2 days can be produced in typical lensing 
systems. A broad range of methods for reliably determin- 
ing time delays from typical data and a deep understand- 
ing of the systematic problems has been developed in last 
years. Time delay estimates are more and more accurate 
and an accuracy of 0.2 days has been a lready obtained in 
the case of B0218+357 jBiggs et alll999l) . so that the detec- 
tion of the gravito-magnetic time delay will be soon within 
astronomers' reach. Observations with radio interferometers 
or HST can measure the relative positions of the images and 
lenses to accuracies ~ 0.005". Shifts in the image positions 
due to a gravito-magnetic field are usually well above this 
limit. 

Our results are nearly unaffected by the presence of 
low-mass satellites and stars. These substructures do not 



have any impact on time delays and can only produce ran- 
dom perturbations of approximate ly 0.001" in image posi- 
tions jKochanek fc Schechteil2004ft . quite below the gravito- 
magnetic effect. Deviations from circular symmetry due to 
either the ellipticity of the deflector or the local tidal gravity 
field from nearby objects should be considered too. However, 
for a singular isothermal model with arbitrary structure, the 
time delay s turn out of be independ ent of the angular struc- 
ture (Kochanc k fc Schechter] 120041) . Other higher-order ef- 
fects, such as the delay due to the quadrupole moment of 
the deflector, should be considered in addition to the gravito- 
magnetic time delay. Unlike other effects, a gravito-magnetic 
field can break the circular symmetry of the lens, ind ucing 
characteristic features in lensing events <Serendl2003bh . At 
least in principle and for some configurations of the images, 
a suitable combinations of the observable quantities can 
be used to remove additional effec ts, due to a quadrupole 
moment, from obse rvational data (|Ciufolini fc Ricci 120021 : 
ICiufolini et~al]|2003h . Furthermore, when the lensed images 
lie on opposite sides of the lens galaxy, the time delay be- 
comes ne arly insensitive to the quadru pole structure of the 
deflector llKochanek fc Schechterll20o4l ■ 

Measurements of gravito-magnetic time delays could of- 
fer an interesting prospective to address the "angular mo- 
mentum problem". Cold dark matter models of universe 
with a substantial cosmological constant appear to fit large 
scale structure observations well, but some areas of possible 
disagreement between theory and observations still persist. 
The most serious small scale pro blem regards th e origin and 
angular momentum in galaxies iPrimackl l2004f) . Two "an- 
gular momentum problems" prevent the for mation of realis - 
tic spiral galaxies in numerical simulations iPrimackll2004f) : 
i) overcooling in merging satellites with too much transfer 
of angular momentum to the dark halo and ii) the wrong 
distribution of specific angular momentum in halos, if the 
baryonic material has the same angular momentum distri- 
bution as dark matter halo. Detection of gravito-magnetic 
effects in gravitational lensing systems due to the spin of 
the deflector, either through time-delays measurements, as 
discussed in this paper, or through observations of the ro- 
tation of the plane of polarizati on of light waves from the 
background source dSerendl200l . could provide direct esti- 
mates of angular momentum and help in developing a better 
understanding of astrophysics in galaxies. 
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